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Abstract 

In recent literature there are an increasing number of papers 
where the forbidden sets of difference equations are computed. We 
review and complete different attempts to describe the forbidden set 
and propose new perspectives for further research and a list of open 
problems in this field. 


1 Introduction 

The study of difference equations (DEs) is an interesting and useful branch 
of discrete dynamical systems due to their variety of behaviors and ability 
to modelling phenomena of applied sciences (see [m Ea EH EH] and the 
references therein). The standard framework for this study is to consider 
iteration functions and sets of initial conditions in such a way that the values 
of the iterates belong to the domain of dehnition of the iteration function 
and therefore the solutions are always well defined. For example, in rational 
difference equations (RDEs) a common hypothesis is to consider positive co¬ 
efficients and initial conditions also positive, see [milHlEHI. 

Such kind of restrictions are also motivated by the use of DE as applied mod¬ 
els, where negative initial conditions and/or parameters are usually mean¬ 
ingless, EH. 

But there is a recent interest to extend the known results to a new frame¬ 
work where initial conditions can be taken to be arbitrary real numbers and 
no restrictions are imposed to iteration functions. It is in this setting where 

^This paper will be published in the special ECIT’14 issue of Journal of Difference 
Equations and Applications 
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appears the forbidden set of a DE, the set of initial conditions for which after 
a hnite number of iterates we reach a value outside the domain of definition 
of the iteration function. Indeed, the central problem of the theory of DEs 
is reformulated in the following way: 

Given a DE, to determine the good (Q) and forbidden (D) sets of initial 
conditions. For points in the good set, to describe the dynamical properties of 
the solutions generated by them: boundedness, stability, periodicity, asymp¬ 
totic behavior, etc. 

In this paper we are interested in the hrst part of the former problem: how 
can the forbidden set of a given DE of order k be determined? In the previous 
literature to describe such sets, when it is achieved, is usually interpreted as 
to be able to write a general term of a sequence of hypersurfaces in But 
those cases are precisely the corresponding to DE where it is also possible to 
give a general term for the solutions. Unfortunately there are a little number 
of DEs with explicitly defined solutions. Hence we claim that new qualitative 
perspectives must be assumed to deal with the above problem. Therefore, 
the goals of this paper are the following: to organize several techniques used 
in the literature for the explicit determination of the forbidden set, revealing 
their resemblance in some cases and giving some hints about how they can 
be generalized. Thus we get a long list of DEs with known forbidden set that 
can be used as a frame to deal with the more ambitious problem of describe 
the forbidden set of a general DE. We review and introduce some methods 
to work also in that case. And finally we propose some future directions of 
research. 

The paper is organized as follows: after some preliminaries, we review 
the Riccati DE, which is one of the few examples of DE where the former 
explicit description is possible. As far as we know, almost all the literature 
where the forbidden set is described using a general term includes some kind 
of semiconjugacy with a Riccati DE. DEs obtained via a change of variables 
or topological semiconjugacy are the topic of the rest of the section. In the 
following we will discuss how algebraic invariants can be used to transform 
a given equation into a Riccati or linear one depending upon a parameter, 
and therefore determining its forbidden set. 

After we will deal with an example, found in [16], of description where 
the elements of the forbidden set are given recurrently but in explicit form. 

We introduce a symbolic description of complex and real points of D in 
section]^ whereupon in section]^ we study some additional ways to deal with 
the forbidden set without an explicit formula. We hnalize with a list of open 
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problems and conjectures. 

To avoid an overly exposition, we have omitted some topics as the study 
of systems of difference equations (SDEs), (see P UHl [HI US [2ll EHl HU Il2l 
ElllS]), the use of normal forms (see |35]), the systematic study of forbidden 
sets in the globally periodic case and the importance of forbidden sets in 
Lorenz maps. 


2 The forbidden set problem 

In this paper we deal with difference equations (DEs) and systems of 
difference equations (SDEs). General dehnitions of these concepts are the 
following. 

Let X be a nonempty set and ^4 C X a nonempty subset of X. Let / : A —)• X 
be a map. A DE of order 1 is the formal expression 


Xn+l = fiXn) ( 1 ) 

([^ represents a set of hnite or inhnite sequences of A, as given xq G A, 
{xn)n=o is constructed by recurrence using ([^, that is, xi = /(xq) and if 
Xn ^ A then Xn+i = f{xn)- When this process can be repeated indehnitely 
then M = +cx) and {xn)n=o i® named the solution of ([^ generated by the 
initial condition xq. 

We call = {xo G A : dn > 0 | /(x„) ^ A or /(x„) is not dehned} the 
forbidden set (FS) of ([^, and Q = A\J^ the good set (GS) of ([^. 

A DE of order k is 

Xn+l f {Xni ■ ■ ■ 1 Xii—k+l) (2) 

where now / : A C X^ —)• X is the iteration map, and the FS is dehned as 

7 = {(Xo, . . . ,X_fc+i) : 3n > 0 I (/(Xn, . . . ,Xn-fc+l),Xn, . . . ,X„_fc) ^ A 

or f{Xn,...,Xn _fc+i) is not dehned} 


In a similar way, we dehne a SDE of order 1 as 


Xn+l ~ 

h{xl,.. 


rpT _ 

■^n+1 ~ 

frixi, . . 


A C X'’ -)■ X, i = 1, . 

. .,r. 



( 3 ) 


System ([^ can be expressed as DE of order 1 of type ([^ using the vectorial 
notation = (x)j,..., xjj) and considering the map F : A C X'' —)■ X'’ whose 
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components are fi,..., fr- 

Finally, a SDE of order k is defined as a set of equations 

\ 

•^n+l Ji\'^ni ■ ■ ■ 1 •^n—k+li ■ ■ ■ i •^ni ■ ■ ■ i '^n—k+lJ 

using maps fi'.Aci X'’*' —)■ X for i = 1,... ,r. 

In vectorial form the set of equations (|^ can be rewritten as in ([^ 

Xn+l = P{Xn, ■ ■ ■ , Xn -k+l) 


with F : ^ XF 


The former definitions depend on how is given the domain of definition 
A of the iteration map. To remark that point, let’s consider c G M \ {0} and 
the real DE 

Xn^l C 

which forbidden set is = A when A is defined as M \ {c} meanwhile F is 
empty when A is the natural domain of the constant function. 


To avoid degenerated or trival cases, some further restrictions must be 
imposed to A. Natural and common restrictions consist usually in regarding 
which is the domain of definition of the iteration map. Therefore, in equation 

1 

Xn+l = - (5) 

Xn 

we say that F = { 0 } as every solution with xq 7 ^ 0 is 2 -periodic and xi is 
not dehned when xq = 0. It is implicitly assumed that A = M \ {0}. 

We remark that ([^ is also a DE over X = MP, the projective line, and then 
F = ^ being every solution 2-periodic if A = X using the rules ^ = cxd, = 0. 


Analogously, the DE 


^n +1 


1 

SfTI 


( 6 ) 


has = 0 in M, i G -F when ([^ is taken over the complex field and F is 
again empty in the projective plane CP. 


In practical applications an undefined zero division means that the 
denominator belongs to certain neighbourhood of zero. Therefore in DE 
as 


^n +1 


P{Xn) 


(7) 
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where P, Q are real polynomials, the forbidden set problem could be studied 
using A = M \ {—e,e) for a machine-value e > 0 . 

In applied models where only positive values of the variables have prac¬ 
tical meaning, A could be 

As a hnal example of different ways to consider the A set, let’s recall the 
DE associated to a Lorenz map. Let / C M be the unit interval [0,1]. Let 
c G (0,1). We say that f : I \ {c} —)■ / is an expansive Lorenz map if the 
following conditions are satished 

i) / is increasing and continuous in each interval [ 0 , c) and (c, 1 ] 

ii) lim f{x) = l and lim f{x) = 0 

x—^c~ 

hi) / is topologically expansive, i.e., there exists £ > 0 such that for any two 
solutions {Xn)n =0 and {yn)n =0 of 0 not containing the point c, there is 
some i > 0 with \xi — yi\ > e, and if some of the former solutions contains 
c the inequality remains valid taking /(c) = 1 or /(c) = 0 . 

Condition (iii) is equivalent to say that the preimages of the point c are dense 
in /, or, in the forbidden set notation, to say that P = /. 

It is obvious that / can be arbitrarily defined in c (without bilateral conti¬ 
nuity) and in that case P = 0. In the standard definition it is assumed that 
A = J\{c}. 

Lorenz maps are an important tool in the study of the Lorenz differential 
equations and the Lorenz attractor, and also in the computation of topologi¬ 
cal entropy in real discontinuous maps. See [^1^ and the references therein. 

In the following we deal mostly with DEs where iteration functions are 
quotients of polynomials, known as rational difference equations (RDEs). 


3 Semiconjugacies on Riccati equations 

Let’s briefly recall some well known results about Riccati DEs. Let a, b,c,d E 
M such that |c| -|- \d\ 7 ^ 0. A Riccati DE of order 1 is 


Xn+l 


a + bxr. 
C -1- dXr 


( 8 ) 


Some special cases occur for particular values of the parameters. If d = 0, 
ad — cb = 0 or 6 - 1 -c = 0 the equation is linear, constant or globally of period 2 
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respectively (see [IHIEH]). If none of those conditions stand, an affine change 
of variables transforms ([^ into 


Un+l — 1 


R 

Vn 


(9) 


R = being the Riccati number. Finally i/n = leads to the linear 

DE Zn+i = Zn — Rzn- 1 - The closed form solution of the linear DE ables to 
write a closed form solution for 0. Moreover, there is a correspondence 
between the solutions of the linear DE containing the zero element and the 
hnite or forbidden set generated solutions of ([^. This is the idea behind the 
characterization of the FS of (|^ that can be found, for example, in [SlfT^EB]. 

From a topological point of view, that characterization shows that is a 
convergent sequence or a dense set in M or a hnite set. Moreover, in the last 
case equation ([^ is globally periodic (see [M]). 


Analogously, for the second order Riccati DE 

b c 

Xn+i = a^ -1-a, 6 , c G M, c 7 ^ 0 (10) 

we can transform it into a linear one using the change of variables . 

The forbidden set is then described as (see [HI [H l39] ) 


+00 

‘T= U {{u,v) \ PinUV +/32nU +Psn = 0 } (11) 

n=—l 


where coefficients {/3in, /d 2 n, /33n}n=-i depend on the roots of the characteristic 
equation. Roughly speaking is a countable union of plane hyperbolas 
convergent to a limit curve in some cases, dense in open sets in other cases 
and even a hnite collection when the DE is globally periodic. 

Our hrst proposal (open problem is to clarify which kind of topological 
objects can be obtained in the Riccati DE of second order and to generalize 
to the Riccati DE of order k 


(2i Go Gu 

Xn+l = ^0 H-^-h . . . H- 

Xn X^iXfi—X Xji . . . 


ao,... ,afc e M, Ofc 7 ^ 0 (12) 


We claim also (conjecture that DE (12) is globally periodic if and only if 
is a hnite collection of hypersurfaces in R^, generalizating what happens 
in the case of order 1 (see [HI]). 

To study the higher order Riccati DE see also [7]. 
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In the former discussion, it was important to determine the set of initial 
conditions whose solutions in the linear equation included the zero element. 
This is called the zero set of the linear equation, Z. For example, in the 
Riccati DE of order 1, Z is defined as the set of points (zq, Z-i) G such 
that the solution of the linear DE includes the zero element. But zq and Z-i 
depend on the initial condition |/o of (|9) via the change of variables yn = 

We can write the dependence as zq = zo^yo), Z-i = Z-i{yo). Therefore, 
can be expressed as 

7 = {|/o e M : 3n > 0 I Zn{y^) = 0} 

= {l/o £ 1^ : {zo{yQ),Z-i{iif))) G Z} 


Generalizing to the Riccati DE of order k (12), we get 


7 = {(xo,... ,a;_A:+i) G 

= {(xo,... ,x_fc+i) G 


3n > 0 I Zn{xQ ,..., x-k+i) = 0} 
(xq, • • • , Z_k) € 2} 


where = Zn{xo ,..., X-k+i) are the terms of the associated linear equation 
of order k + 1. 

Note that the semiconjugacy formula yn = has zero as a pole. The 
forbidden can be regarded as the transformation of certain special set of the 
linear equation corresponding to the singularities of the semiconjugacy. 


There is an increasing number of works in the recent literature that use 
the closed form solution of the Riccati equation to describe the forbidden 
set of other rational difference equations. Let’s review and complete some of 
them. 


The order 3 DE 


^n +1 — 


(XXyiXYl—\ 

hXn + CXn-2 


(13) 


semicon- 


can be transformed into the linear form Zn+i = - + -Zn using the 
jugacy Zn = The forbidden set is a sequence of planes in The case 
6 = 1 and a,c > 0 in (13) was developed in [38], while a = ±1, 6 = 1 and 
c = —1 is in [1]. See also mg for a > 0 and he < 0. 

A similar change of variables, Zn = gives the former linear equation 
when it is applied to the following DE of order 4, see [2]. 

ClXfiXYi—2 


^n +1 — 


a, 6, c > 0 


(14) 


hXn + CXn-3 

We propose a generalization of this problem remarking that the following DE 

(15) 


aXnXn—k+l / ^ 1,1 1 1 ^ p. 

Xn +1 = ^-;- a [J, \b\ + |c| > (J 


hXn T CXn—k 
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becomes Zn+i 


= a + with Zn = ^ (open problem |^. 


In [ 13 ], the change Zn = XnXn-i gives the Riccati eqnation Zn+i = 
when we apply it to 

Xn—l 

^n+1 I 7 

a + bXnXn-1 

The anthors explicitly describe as a seqnence of plane hyperbolas. 


Zn 

a+bzn 


(16) 


Eqnation (16) admits several generalizations. If we nse the change of 


variables Zn = XnXn-k applied to the Riccati DE Zn+i = 


a-\-bzn 


we obtain 


^n^n—k 

^n +1 I 

(lXn—k-\-l ^n^n—k-\-l^n—k 


( 17 ) 


a RDE of order A; + 1 whose forbidden set can be described using that of the 
Riccati DE. In this case we have a RDE of third degree. A more interesting 
family of equations, where the degree remains equal to two, is 






Qj “ 1 “ hxYi—iXji—2i—l 


(18) 


where i is a natnral nnmber. Remark that for i = 0 we get eqnation (16). 
Using the change of variables Zn = Xn-i-iXn, (18) leads to 

^n—i 


^n +1 


a + bzn-i 


that it is not a Riccati difference equation, but can still be reduced to a linear 
form. Indeed, if with an affine change we transform it into 

-1 ^ 

Vn+l — 1 

Vn—i 


we can introduce ?/„ = — to get 

'R'n —2 —1 

^n+l '^n—i 2i—1 

a linear DE of order 2i + 2 (see open problem]^. 


In [30], it is shown that R is a family of straight lines in in the case 
of equation 

aXn-l + bXr, 


CXn-1 + dx, 


-Xr, 


^n +1 — 


(19) 















via the change of variables 


Here we get the Riccati form Zn+i = 

y _ Xn 

1 

Given i > 0 a possible generalization is 


a-\-hzn 

C+dZn 


^n+1 


QjXn—2i—\ ^^n—i 

CXn—2i—l ^^n—i 


^71—i 


( 20 ) 


Since Zn = 


^n — i — 1 


rednce (20) to 


Zn+l — 


a + bzn-i 

C ~\~ djZji—i 


that can be linearized to an eqnation of order 2z + 2 as we do for 
problem . 


(open 


The difference eqnation 


^n +1 


^n—2 

itl “t“ XjiXfi—\X'ii—2 


( 21 ) 


of reference BTl becomes Un+i = with the change Un = XnXn-iXn- 2 - 

Generalizing, in [10] the RDE of order k + 1 


^n +1 


(y.Xn—k 

13 + JXn . . . Xn-k ’ 


a, /S, 7 > 0 


( 22 ) 


is rednced to the Riccati form 


^n+l 


aZn 


f3 + 'yZr, 


(23) 


when we nse the mnltiplicative change of variables Zn = Xn ■ ■ ■ Xn-k- 
Let’s generalize the FS characterization of m to the case where a, (3,'y 


are arbitrary complex nnmbers. If one of them is zero, then eqnation (22) 


becomes trivial or globally periodic and 3F is empty or the set of points with 
at least one zero component. 


If af3'y 7 ^ 0, let c = §. Therefore (23) transforms into 


^n+l 


CZr, 


1 + Zn 


The former eqnation admits the following closed form solntion 


^n. — 




1 + (1 + c + ... + c^~^)zo 
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And from here it is easy to give the forbidden set expression. If c = 1, then 
Zn = and Ifcis another root of the unity, then the 

equation is globally periodic and the forbidden set of the Riccati equation is 
hnite: 
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1-c 1-c ) 

7 — 1 ’ ’ ~ 1 J 


supposing m to be the smallest positive integer such that c”* = 1. 
When we apply these results to (22), we get that 
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M ^ 

U < (x-fe,..., Xo) e : X-fc ... Xo 

n=0 



(24) 


that is, a countable union of generalized hyperbolas when c is not a root of 
the unity (M = +cx)) or a finite union if = 1. 


There are a number of RDEs for which the closed form solution and FS 
is given in [32]. Those equations can be grouped in three categories: 

• Those obtained from the multiplicative DE 


Xn+i = ax^xl_^, a 7 ^ 0 p,q eZ (25) 

when we do change of variables as 


Vn b 


Vn + b 

Vn + C 


We deal then with RDE of order 2 and 
to obtain the closed form of (25) from 
constructed. 


degree 2. It is not difficult 
which the FS expression is 


• RDE of order 2 and degree 2 resulting of the introduction of the variable 
ffi the general linear DE 


Xn-\-l dXn T bXn—\ T C 


(26) 


• And another family of RDE of order 2 and degree 2, given in this case 
by the change applied to the Riccati DE of order 1 


aXn + b 

Xn+i — I y 

CXn + a 


(27) 
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In each of the former cases we propose some generalizations. Given 
a Mobius transformation T[x) = it must be possible to explicitly 

determine the FS of the RDE constructed with the change = T(|/„) applied 
to equation (25) or applied to equation (26). These are the claims of open 
problems and respectively. 

In the third case, we have that the closed form solution of (27) and the 
change of variables able to write 


Xn = g{a, b, c, d, yo, y-i) = gn 

Moreover, the relationship Xn = can be interpreted as a 

nonautonomous linear DE when solved for j/n+i, that is 


Vn+l — 


^g-n ^ f^gn (d 

- yn+ - 

1 - 75'n 1 - 7fi'n 


Therefore a general expression for y^ and for can be computed, 
and of course, the same idea must work for every change of variables 
Xn = Hijjn+ii ■ ■ ■ iyn-k) where the explicit solution of the nonautonomous 
yn equation is known, and for every equation Xn+i = G(xn, ■ ■ ■, Xn-k) whose 
closed form is also known (open problem . 


4 Use of invariants 


An interesting modihcation of the former ideas is the use of invariants to 
describe the FS of some RDEs. Consider the following example from 


Xn+l 


Xr, 


1 + BXn-l - BXn 
where B E C \ {0}. This equation has the following invariant 


-h S ) (1 + BXn-l) — C 

^T). 


(28) 


(29) 


that is, for every solution of (28) there exist C = G(a;o,x_i) such 


that (29) holds for every n > 0. 


The presence of an invariant ables to write an alternative form of the DE. 


Indeed, solving (29) for Xn and changing the indices, we get the following 
Riccati DE 

1 + BXn 

^ C-B- B^Xn 
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where recall that constant C depend on the valnes of x_i and xq. This is the 
basic idea to describe the FS of (28). 

Note that formnla (29) implies the identity 

+ i? ) (1 + BXn) = (-h i? ) (1 + BXn-l) 


^n+1 


Xr 


from where (28) is dednced. 

A generalization of the former remark is the following. Let Ti and T 2 be two 
Mobins transformations, Ti{x) = * = 1,2, and consider the invariant 


Therefore 

and from here we get 


Ti{Xn)T2{Xn-l) = C 

Ti{Xn+l)T2{Xn) = Ti{Xn)T2{Xn-l) 


(30) 


Xn+1 


f Ti{Xn)T2{Xn-l) 


(31) 


T2{Xn) 

a RDE of order 2 whose FS conld be described with the former methodology 
as the invariant (30) implies also that 


Xn+l — 


-1 


C 


T2{Xn) 


which is a Riccati DE of order 1 (open problem]^. 

For example, given Bi and B 2 non zero complex nnmbers, the invariant 
{xn + Bi){xn-i + B 2 ) = C leads to DE 


Xn-\-l 


XnXn-l + {B 2 - Bi)Xn + BiXn-1 


Xn + B2 


and the invariant = C to DE 

Xn-l + l 


Xn-\-l — 


xl + (1 + Bi)Xn - BiXn-1 
Xn-1 + 1 


both of them not inclnded in reference ED. 
Also, given k > 1, the invariant 

Ti{Xn)T2{Xn-k) = C 


(32) 
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produces the DE 


^n+l 


— 1 f T\{x^T2{Xn—k+l) 


T2{,Xn) 

that can be studied in the same way because the C dependent DE 

C 


(33) 


Xn+l — 


-1 


^^2(^71 —fc + lj 

is reduced to a linear equation of order k + 1 using the same changes as in 
equation (18). 


Another example of use of invariants is in |S]. DE 

Xl + Xl_^ - Xn{Xn-l + Xn-2) 


^n+1 


^n—1 ^n—2 


(34) 


verifyes that /(a:„,x„_i) = ^ is constant over each solution of 


^n—1 

(34). Therefore, from the equality 

Xn+l (^n—1 

Xn+1 Xyi 

the following linear relation is deduced 


= C 


(35) 


Xn+1 


c 


-Xn - 


C-I " C -I 


Xn—1 


(36) 


and from here it is easy to give the closed and forbidden set of (34). In 
particular, the FS is 

F = {{x^y^ z) ■. X = y oi y = z} 

A possible generalization is to consider invariants of the form 

^n+1 ^n—k 


= C k,l€N, I 


(37) 


^n+1 ^n—l 

as every solution will be associated to a linear recurrence (open problem]^. 


Finally, a general question is to determine which relationship exists 
between the poles and zeros of an algebraic invariant, from one side, and 
the elements of its associated DE, from other. 
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5 Forbidden set curves given in an explicit 
form 


One of the oldest examples in the literature concerning the forbidden set 
problem is in [16]. Let p < — 1 be a real number, and consider the following 
RDE of order 2 

= p + ^ (38) 


Map F{x,y) = ^p,p+^j is the unfolding of (38) as we have that 

F{xn-i,Xn) = {xn^Xn+i), V?T. > 0. Let G{x,y) = {x{y — p),x) be its inverse 
map, and A = {(a:, 0) : a: G M}. Therefore 


+ 00 


= U G“(A) 


n=0 


This is an obvious characterization of the FS as the set of inverse orbits of 
poles of the iteration map that can be improved as follows. Let’s consider the 
subset y4+ = {(x, 0) : X > 0} and let {gn)t=i be the sequence of functions : 
[0, +CX)) —)• M defined inductively using the operator hg{x) = x- (—p + p“^(x)) 
in the following way 


Pi(x) = -px, gn+i{x) = hg^{x) 


Therefore we have 


Theorem 1 ([T6| theorem 2). Let he the subset of the forbidden set of 
defined by F+ = UnroG'"(R+). Then 

= A^U {(x, y) e [0, +oo) xR:3n>l,y = gfi^{x)} (39) 


In general, given a RDE we can describe its FS as a set of implicitly 
defined hypersurfaces. In the second order case we get a set of implicitly 
defined curves. For example, in the case of Pielou’s equation 

aXn 

Xn+1 = —- 

1 + Xn- 



we can consider the unfolding map F(x, y) = (p, such that 

F(x„_i,x„) = (x„,x„+i). By iterating map F we get 


F^{x,y) 


ay 


afy 


1 + x’ (1 + x)(l + y) 
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F^{x,y) = 


a^y 


a^y 


F^{x,y) = 


a^y 


(1 + x)(l + ?/)’ (1 + |/)(1 + x + a 2 /) 

a^{l + x)y 


(1 + ?/) (1 + X + ay) ’ (1 + X + ay) {1 + x + y + a'^y + xy) 


From where the following hrst forbidden curves are deduced (see hgure ; 

F = {(x,y) e : 1 + X = 0} U 

{(x,y) e : 1 + ?/ = 0} U 

{(x, y) e : 1 + X + aj/ = 0} U 

{{x,y) e : 1 + X + y + a?y + x?/ = 0} U 


Remark that those curves are given in an implicit form. We propose to find 
a way to construct the explicit form as in theorem]^ (open problem [Io| . 


6 Symbolic description of the FS 

Let’s consider the following RDEs of order 1 and degree 2 


^n+1 



(41) 


^n+1 



( 42 ) 


where a is a real or complex parameter. These equations are studied in [12] 
in connection with their relationship with Li-Yorke chaos in RDEs. 

It is possible to describe the forbidden set of the particular case 





using a symbolic notation. Let 


h+(x) 

h-{x) 



( 43 ) 


( 44 ) 


be the inverse branches of the iteration function in (43). 
forward multiorbit of the poles of /(x) = Moreover 


Clearly, F is the 
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Figure 1: The first 10 forbidden curves in Pielou’s RDE. From up to bottom 
and left to right, a = 1, a = 2, a = —1 and a = —2 


Proposition 6.1. The unique preimage of pole —1 in (fS) is 0. The 
preimages of pole 1 are ai... where Oj = or a* = h._ for i = 1,..., n, 
and ai.. .On is the abbreviation of ai o ... o a„(l). Those preimages are 
always well defined in C, that is, Oi... a„ ^ 0 for every word, and the 
former representation is unique, that is, given two words such that oi... a„ = 
bi.. .bm then n = m and a* = bi for i = 1,... ,n. 

Therefore 


T = { — 1, 0,1} U {oi ... On ■ n E N,ai = h±} C C 


(45) 


Proof: Let’s see that words ai... a„ are distint from zero and are unique, as 
the remaining assertions are straightforward to demonstrate. 
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When n = 1, ai(l) = ±a/2 7 ^ 0. For n > 1 , if oi... a„ = 0 , then 
ai(a 2 ...a„) = 0 and therefore 02 ... a„ = — 1 , what is impossible as — 1 
does not belong to the image of h_|_ or h_. 

Now, suppose that oi... = 61 ... bm- There is no lack of generality if we 

assume that oi 7 ^ &i, because and h- are injective maps. Let x = 02 ... 
and y = 62 • • • &m- The former equality is rewritten as ai(x) = bi{y). But 
note that Imh+ fl Imh_ = { 0 } from where we deduce that x = y = —1, 
what is again impossible. To avoid the contradictions both words ai... a„ 
and bi.. .bm must be equal. □ 


Note that this is a set in the complex field (see figure [^. If we center 
in the real field, note that applying h- one time we stay over the reals, and 
applying it two consecutive times we get C \ M. Therefore 


Corollary 6.2. 


The forbidden set of (43) over M is 


J'R = {“Ij 0 , l}U{ai... a„ : n G N, Oj = h±, there is no two consecutive h_} 


Numerically, Jtr is the projection of figure over the real line. In open 



Figure 2: Forbidden set of x^+i = 

J- 


problem we propose to find the symbolic description of 4Fc and Jtr for the 
families of RDEs (41) and (42). 
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7 Qualitative approaches 

As we have seen in the former sections, there are some works describing the 
FSs of DEs. But it should be remarked that the most part of DEs don’t have 
a closed form and therefore there is probably no explicit expression for their 
FSs. 

We claim that new perspectives should be taken in order to solve the problem 
of the FS for a generic DE. Specihcally, given a DE we propose three general 
questions about its FS 

1. Which are its topological and metrical properties? (Openness, 
Boundedness, Fractal dimension,...) 

2. What can be said about perturbations of its parameters? 

For example, given a family of DEs Xn+i = fa{xn) depending on 
parameter a, and supposing that some feature of the FS is known 
when a = Oq, what can be said about the FS of the DE when 
a, G (no — A 

3. Are there some relationships between the properties of the FS and the 
dynamical behavior of the DE? 

There are, as far as we know, very few papers with a qualitative approach. 
Let’s briefly review them. 


In [M] the problem of closedness of the FS is discussed in connection 
with the global periodicity question. Let IK be the real or complex held. Let 
F : —)■ and consider the fc-dimensional DE or dynamical system 


Xn+l = F{Xn) 


(46) 


Let D C IK*’ be the domain of F. A set A <Z D has full Lebesgue measure in 
D if the measure of D\ A is zero. The map F is almost a local diffeomorphism 
if there exists an open set V C IK*’ of full Lebesgue measure in D such that 
F is continuously differentiable on V (or analytic on V when IK = C) and 
DF{x) 7 ^ 0 a.e. x eV. 


Theorem 2 ([31], theorem 3.). Let F : D C IK*’ —)■ IK*’ 6e an almost local 
diffeomorphism such that the set F“*(IK*’ \ D) has zero Lebesgue measure. 
Suppose also that D is open and that F is a continuous globally periodic 
map. Then the good set of (46) is closed and has full measure in D. 


There are some consequences of this result. The natural domain of a 
rational map F = (Fi,...,Ffc) is the set of points in IK* such that any 
denominator in F vanish. The hrst result is for the rational case. 
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Corollary 7.1 

map F : D 




, corollary 4). Let D be the natural domain of a rational 
—)■ If F is globally periodic and detDF{x) is not 


identically zero, then the good set of (46) is open and has full Lebesgue 
measure in 


In the case of RDE, we say that a function f{xi,... ,Xk) depends 
effectively on its last variable if 7 ^ 0 for some x = (xi,... ,Xk) G D. 

Therefore for RDE 

Xn+l = f{Xn, Xn-k+l) (47) 

we get 


Corollary 7.2 ([31], corollary 5). Let f : D C ^ K be a rational 
function, and D its natural domain. If f depends effectively on its last 
variable and (fl) is globally periodic, then the good set is open and has full 
Lebesgue measure in K*'. 


A more concrete application is given in the case of the Riccati DE of order 1 

Proposition 7.3 ([31], proposition 7.). DE ^ such that none of the 
quantities d, ad — cb and b + c is zero, has an open good set if, and only 
if, the equation is globally periodic. 


It is a natural question to wonder if this result applies also to any RDE. 
Again in [31| is given the following not globally periodic example 


^n+1 


{Xn l)(^n 2) 

Xn 


(48) 


whose FS is = {0,1, 2} and therefore closed. Numerical experiments show 
also that F is not closed when the equation is regarded in C (see figure 
open problem [I^ and conjecture [^. 


A second example of qualitative study of DE is in [33] • There the 
author stands sufficient conditions for a DE to have a good set of full 
Lebesgue measure, generalizing a previous remark of [201 section 2]. In the 
nonautonomous case, the main result is 


Theorem 3 ([33|, theorem 5). Let f : D F K be a function locally in 
fact of order k and such that the set M = {x E D : {f{x),X 2 , ■ ■ ■ ,Xk-i)} has 
Lebesgue measure zero. Then the good set of equation (fl) has full measure 
in D. 
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Figure 3: 50000 elements of the FS of (48). Green points are the unique 
elements of in M. Fixed points are marked in yellow, and they seem to be 
in the clousure of T. 


Here a function / is locally in fact of order k if there exists an open set 
V C D of full Lebesgue measure in D and such that 7 ^ 0 a. e. x eV . 

In [33] it is remarked the importance of determining the topological and geo¬ 
metrical properties of the good set in terms of the iteration function and its 


domain, and two related open problems are proposed (see open problems 13 
and[T4j). 


Also in |l3] there is a generalization of these considerations, studying the 
problem of avoiding forbidden sequences. 


In [3n] the following DE are considered: 


a 

Xn+l — + 1 

(49) 



a 1 

(60) 

^n+1 q p 


■^n-1 
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where a ^ 0 and p and q are odd rationals with at least one of them positive. 
An odd rational is a fraction of the form |^, i,j G N- 

Again the objective is to determine somehow the forbidden set. The 
anthor uses the term crash set referring to T, emphasizing that the discon¬ 
tinuity is a pole. 


None of the forbidden or crash sets of the former equations is explictly 


computed, but it is described qualitatively. In the case of (49), the iteration 


function belongs to the family of monotonic maps with pole, dehned as the 
class of functions / : M —)■ M verifying 

(Al) / is continuous on its domain M \ {0} 

(A2) / is injective on M \ {0} 


(A3) 0e/(M\{0}) 

(A4) lini |/(a;)| = cx> 

Note that / must be an increasing or a decreasing map, and that consists 
in the collection of preimages of the pole 0 , qualitatively described in the 
next theorem. P stands for the set of hxed points of /. 

Theorem 4 ([36], theorem 1). Let f be a monotonic map with pole. 

(a) If f is a decreasing map, let x be the fixed point in (/“^(O), 0). Then the 
sequence T = (/“"'(0))„ is contained in the interval [/“HO); /~^(0)]; 
converges either to x or to a 2-cycle. 

(b) Suppose that f is an increasing map, P 7 ^ 0, and x = inf P > /“^(O). 
Then P is contained in [/“^(0),a;) and converges monotonically to x. 


Remark that there are some cases not included in this result. When / 
is increasing and x = supP < /“^(O), cobweb diagrams show that P is a 
subset of {x, /“Ho)] ^^0 converges monotonically to x. Whereas for P = ^, 
P could be hnite or could have a more complicated structure. For example 
in [36] is wondered if in the last case P can be dense outside a compact set. 
(See open problem [Ts] and £gure|^. 


Theorem 1^ applies not only to equation (49), but also to equations as: 


Xn+l — - 1 " 1 

smh Xn 


( 51 ) 
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Figure 4: Forbidden sets of function /(x) = 1 — ^ represented with a cobweb 
diagram of the inverse iteration starting at the pole 0. From left to right and 
up to bottom, p equals to 3, | and | respectively. 250 iterations were 

taken in each case. 


and, in general to DEs of the form 

^n+l 

where a 7 ^ 0 and 0 : M —)■ 


+ 1 


is a bijection such that 0 ( 0 ) = 0 . 


( 62 ) 


The FS of (50) is given in [36] in terms of the values of a and p. We 
propose to make a similar study for the family of DE 


Xn+i=T{xPJ (53) 

where p is an odd rational and T is the Mobius transformation T(x) = 
with a, 6 , c 7 ^ 0 (open problem [l7|) . 

In the case of equation ( [50| , in [SB] there is a description of the FS curves 
in an explicit form using a recurrent algorithm, in a similar way as it was 
made in section [5l 
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Additional examples of estimative approaches where the FS is graphically 
represented can be found in [TJl |T5l [19] for the following DEs 


1 

^n+l I 

Xn "h Xji—2 


Xn—1 

Xn+l — Z I 
1 

In a similar way, we can estimate the FS of 

A B 

Xn+l — -1- 

Xn Xn—1 


(54) 

(55) 


(56) 


for differents values of parameters A and B. See hgure and open problem 

m 


8 A list of open problems 


Open Problem 1. 

in the second order 

(H- 


Clarify which kind of topological objects can be obtained 
Riccati DE and generalize to the Riccati DE of order k 


Conjecture 1. DE (12) is globally periodic if and only if D is a finite 
collection of hypersurfaces in 


Open Problem 2. Describe the ES of RDE (15). 


Open Problem 3. Describe the ES of RDEs (11) and (18). 


Open Problem 4. Describe the ES of RDEs (20). 


Open Problem 5. LetT{x) = Describe the ES of the RDEs obtained 

by applying the change = T(yn) to DE 


Open Problem 6. LetT{x) = Describe the ES of the RDEs obtained 

by applying the change Xn = T(yn) to DE 


Open Problem 7. Let Xn+i = G(xn, ■ ■ ■ ,Xn-k) a DE of order k + 1 such that 
its closed form is known. Let Xn = H{yn+i, ■ ■ ■, yn-k) o, change of variables 
that can be rewriten as a nonautonomous DE yn+i = H{yn, ■ ■ ■ ,yn-k) with 
coefficients depending of Xn and such that its closed form solution is also 
known. Determine the forbidden set of the DE obtained by applying the 
change of variables to Xn+i = G{xm • • •, Xn-k)- 
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Open Problem 8. Describe the FSs of DEs (31) and (33) having algebraic 
invariants. 


Open Problem 9. Let k and I be distint natural numbers greater than 0. 
Determine the closed form solution and the forbidden set of DEs admiting 
an invariant of the form [31). 


Open Problem 10. For every RDE of order 2, determine a region A C 
such that the forbidden curves in A can be given in explicit form, and 
determine them. 


Open Problem 11. Find the symbolic description of Dc and for the 


families of RDEs (fl) and (42) 


Open Problem 12 (|3l], open problem 1). Let F : D C ^ an 
almost local diffeomorphism such that F~^{M.^\D) has Lebesgue measure zero. 
Suppose that eguation (46) is globally periodic. Obtain sufficient conditions 
in order that 4F is closed. 


Conjecture 2 ([3l], conjecture 2). Let f : D ^ C be a rational function and 
D its natural domain. Assume that f depends effectively on its last variable. 
Then the good set of (41) is open if, and only if, the equation is globally 
periodic. 


Open Problem 13 (|33], open problem 1.). To obtain necessary conditions 
in order to have that an equilibrium point of an autonomous difference 
equation belongs to the interior of the good set. 


Open Problem 14 ([33]. open problem 2.). To obtain necessary conditions 
in order to have that the good set of a difference equation is open. 


Open Problem 15. Let / : M —)■ M 6e a monotonic map with pole. Complete 
the results of theorem\^ in the sense of describing qualitatively the FS of DE 
Xn+i = f{xn) when f is an increasing map with no fixed point. In particular, 
answer to the following questions: 


1. In which cases IF is finite? 

2. Can be IF dense outside a compact set? (See JEBf)- 


Open Problem 16 (|2T]). Determine the FS of DE (56). 

Open Problem 17. Determine the FS of DE (53) by applying theorem^in 
terms of a, b, c and p. 
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As we have seen throughout the paper, there are a lot of DEs with closed 
form solution and explicit forbidden set. There are many ways to generalize 
them and therefore the literature on this subject is growing with more and 
more works devoted to that kind of DE. Although those works have the 
value to enrich the set of known examples, may be it is time to put some 
order in the field. We propose to elaborate a database of DEs and SDEs 
including, among others, the following topics: the forbidden and good sets, 
the asymptotic behavior, the boundedness of solutions, the relationships with 
other DE or SDE and the applied models related with them. 
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Figure 5: Two representations of the FS of DE (56). On the right the first 
10 forbidden curves are drawn. On the left we color after the value of the 
tenth iterate of each point, reserving gray for points of . From up to down, 
the values of the parameters (A, B) are the following: (1, 5), (1,1), (—1, —1). 
Numerical experiments show also tftMi when {A^B) equals to (1,-1) and 
(—1,1) respectively, vertical symmetries of the second and third row are 
obtained. 


























































